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This paper discusses a method for numerical evaluation of the vibrations 
of a cylindrical shell structure induced by a low-speed external turbulent flow. 
The direction of flow is along the axis of revolution of the shell (see Figure 1), 
and the source of excitation is the pressure fluctuations in the turbulent bound- 
ary layer (TBL). 

For the investigation of vibration and noise problems it is usually more 
desirable to utilize the modal expansion approach. The axisymmetric shell struc- 
ture shown in Figure 1 can be modeled by the assemblage of conical-shell finite- 
elements. This modeling allows the eigenfunction ^(x.O) to be represented 
in a rectangular product of a longitudinal modal function fnmCx) and a circular 
harmonic function cos ra0 (or sin m0), i.e., 

^mn^ x »®^ — fmn( x ) cos ra0 (1) 

m=0 ,1,2,... 

n=l, 2, 3, . . . 

The forcing function from the TBL is assumed to be spatially homogeneous 
and temporally stationary. It is commonly expressed in terms of the wavevector- 
frequency spectrum $p(k 1 ,k 3 ,oj) where k^ is the streamwise wavenumber and k 
is the transverse wavenumber. For the calculation of the structural acceptance 
this forcing function, the structural modes must also be expressed in wave- 
number space. This can be accomplished by taking a spatial Fourier transform 
of the modes. The finite-element modeling provides the computed eigenfunction 
defined at a set of discrete points. If the grid points on the flow surface are 
equally spaced, a Fast Fourier Transform (FFT ) routine may be used. From the 
FFT spectral coefficients, we may express 
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f mn' x ' “ a o + L £a v cos x + b v sin x] (2) 

v=l L L 

where U is one less than one half of the total number of FFT data points and L 
is the axial length of the structure. 

The effective modal input spectral density ^(w) from the TBL can be 
evaluated as follows, 

x mn( u > “ A 2 + <*>) J mn,ran< u ) (3) 
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where A is the total area of the flow surface, Jnmjmn^) * s tbe moc * ab j°i nt “ 
acceptance (or called self - acceptance) , which provides a measure of the degree of 
coupling between the turbulent pressure field and the structure, and <J>(w) is 
the frequency spectrum of the TBL pressure fluctuations. The joint— acceptance 
can be computed by the following summation, i*e., 

J mn,mn( u ) = [ (2 u) 3 /A<J>(o)) ] 


m 1 Jj „ 2vn m .. 

[a§ $ p( 0 »r» u ) + — l ( a v +b v) $ p ( ~ > 7 

r R 2 v== ]_ L K 


where R is the radius of the cylindrical shell. 

Often it is required to evaluate the summation up to the convection wavenumber. 
This requires that the number of FFT data points to be approximately u>L/7rtJ c , 
where U c is the convection velocity of the TBL. If the length L of the structure 
is large and the frequencies of interest are high, the required data points will 
generally exceed the number of finite— element grid points. This difficulty can 
be overcome by obtaining additional data points from spline fitting and interpola- 
tion of the eigenfunctions. 

If the size of structure is larger than the correlation length of the pressure 
field, the cross-modal acceptance Jnm,m T n f > can ne gl ecte d* this case, 

the structural displacement response spectrum [s(m)] evaluated at 0 = 0 O can be 
calculated as follows: 

[ s (to) ] = ['l»][H^(io)JfTT im ((o)J[H mn (to)][t)t] T (5) 

where [iji] is the assembly of eigenvectors, each column represents one eigen- 
vector {f ran ( x i) cos m6 o}’ H mn (to) is the modal admittance function which is 
defined as 

H mn (w) = i M mn[ (a, mn _a)2) + i((0 ran r W + “ a, mn 6 mn>H 1 (6) 

and where M m is the mode mass, is the structural modal loss factor, and 

6 ran is the modal acoustic loss factor. 

The most difficult task in the numerical evaluation of flow induced vibration 
is the uncertainty about the forcing functions, i.e., the wavevector-f requency 
spectrum. Several theoretical forcing function models have been published in 
recent years, all of them require an empirical fit with experimental data. Pub- 
lished experimental data are very widely scattered depending on the measuring 
facility, surface property, the method of scaling the data, etc. Selection of a 
suitable forcing function model and the experimental data thus depend heavily on 
experienced engineering judgement and knowledge of how the experimental data are 
obtained. 
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FIS. 1 A CYLINDRICAL SHELL STRUCTURE WITH 
AN EXTERNAL TURBULENT FLOW 
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